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1.  Introduction.  Many  physical  problems  require  the  solution  of 
partial  differential  equations  on  some  infinite  domain  ft  with  boundary 
3ft.  An  example  is  given  in  fig.  1. 


^ Jn 


ft 


x*U~ 


fig.  1. 


For  computational  reasons  ft  is  replaced  by  a finite  domain  ft 


B, 


v/-1 


x*0 


ft  - ft. 


x=a  B, 


fig.  2.  ^ 

and  the  problem  arises  to  specify  boundary  conditions  at  the  artificial 
boundary  Bg.  Consider,  for  example,  the  differential  equations  for  a 
nonviscous  fluid  which  at  subsonic  speed  leaves  ft^  through  the 
boundary  Bg.  Then  there  is  one  characteristic  which  points  back  into 
the  region  ft^  and  therefore  one  boundary  condition  has  to  be  given. 

In  general,  no  detailed  knowlecfes '*©f  the  solution  on  Bg  is  given  and 
other  principles  have  to  be  applied.  For  example,  if  one  solves  the 
problem  by  a difference  approximation  then  one  often  uses  upstream  diffe* 
rencing  on  Bg  for  all  the  dependent  variables  (see  for  example  Roache 
[k]).  This  procedure  is  sometimes  combined  with  overspecifying  the 
solution  on  the  boundary  B^. 

Recently,  B.  Engquist  and  A.  Majda  [1]  have  proposed  another 
aciple , namely,  to  specify  the  boui 
a way  that  no  reflection  takes  place. 


principle,  namely,  to  specify  the  boundary  conditions  on  Bg  in  such 
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In  this  paper  we  want  to  investigate  when  these  principles 
garantuee  that  the  solution  of  the  simplified  problem  is  close  to 
the  solution  of  the  original  one.  A necessary  condition  for  this  is 
that  in  ft  - ft^  the  solution  of  the  original  problem  only  changes 
slowly  with  respect  to  space  and  time.  Therefore  we  can  linearize  the 
problem  and  we  assume  that  the  linearized  equations  represent  a 
hyperbolic  first  order  system.  This  is  true  for  ideal  flow  problems. 

In  the  next  section  we  consider  the  model  equation 


(1.1) 


3u/3t  ■ 3u/3x  , x ^ 0,  t >_  0, 

u(x,0)  «.f(x)  for  t = 0, 


on  the  half  line 
3v/3t 

(1.2) 

v(x,0) 


0 £ x < » and  approximate  it  by 
' 3v/3x,  0£x£a,  t > 0, 

» f(x)  for  t « 0, 


on  the  finite  interval  0 <_  x £ a.  The  last  problem  is  well  posed  if 
we  specify  boundary  conditions  at  x * a (but  not  at  x * 0). 

We  solve  (1.2)  by  the  Lax-Wendroff  difference  scheme.  As  boundary 
conditions  for  the  difference  approximation  we  have  a number  of  diffe- 
rent possibilities. 

a)  For  x * «e  either  specify  v,  for  example, 
v(0,t)  ■ 0, 

or  we  use  an  extrapolation  procedure 

(hD+)Pv(0,t)  “0  , hD+v(0,t)  ■ v(h,t)  - v(0,t). 

Far  p ■ 2 this  is  the  usual  upstream  differencing. 

• b)  For  x ■ a we  have  the  analogous  possibilities. 

We  obtain  the  following  results. 
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1)  For  x = 0 one  shall  not  overspecify,  i.e.  v shall  not  he 
given.  Upstream  differencing  gives  good  results. 

2)  If  f(x)  * const,  for  x a then 

(1.3)  hD_v(a,t)  « v(a,t)  - v(a-h,t)  « 0 

can  be  used.  Upstream  differencing  at  x * a can  give  completely 
wrong  results . 

3)  The  principle  of  no  reflection  at  x * a is  useful  only  if 
f(x)  * 0 for  x ^ a. 

In  section  3 we  show  that  the  corresponding  conclusions  hold  for 
systems 

3v/3t  * A 3v/3x 

where  A is  a constant  matrix. 

In  section  U we  make  a thorough  investigation  of  systems  with 
variable  coefficients 

3u/3t  * A(x,t)3u/3x. 

We  show  that  extrapolation  procedures  or  the  principle  of  no  reflection 
is  useful  only  if  A(x,t)  is  essentially  independent  of  x.  for  x ^ a. 
If  one  knows  the  asymptotic  behavior  of  A(x,t),  for  example, 

A(x,t)  ■ AQ  + x ^(x.t) 

then  one  can  derive  new  principles,  which  are  useful  for  steady  state 
calculations. 

The  last  section  is  concerned  with  systems 
3u/3t  * A 3u/3x  + B 3u/3y 

in  two  space  dimensions.  The  above  principles  are  useful  if  B 3u/3y 
is  small,  i.e.  the  problem  is  essentially  one  dimensional.  For  the 
ease  that  B is  not  small  we  again  construct  new  principles. 


In  many  applications  the  time  dependent  equations  are  ured  to 
obtain  the  steady  state  solution.  We  shall  also  study  the  effect  of 
our  boundary  conditions  on  the  convergence  rate  to  the  steady  state. 
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2.  The  model  problem.  We  want  to  solve  the  differential  equation 

(1.2)  by  the  Lax-Wendroff  scheme.  Let  k > 0,  h * a/N,  N natural 
number,  denote  the  time  step  and  the  space  step  respectively.  The  grid- 
points  are  given  by  r^  » vh,  tQ  = nk,  v«*0*l,. ...  ,N;  n=0,l,2,...  ; 

and  the  gridfuntions  by  w^n  = v(^,tn).  We  approximate  (1.2)  by 


vvD+1  * (l+kD0+^-D+Djvvn  , 0 < V < N, 


(2.1) 


vv°  ■ fv  • 


where 


2hVv 


v+1 


- w 


v-1’ 


v+1 


~ V 


* w. 


‘ wv-l  * 


denote  the  usual  centered,  forward  and  backward  difference  operators 
respectively.  We  assume  that  k/h  <1,  which  grantuees  that  the  approxi- 
mation is  dissipative. 

We  first  study  the  case  that  proper  boundary  conditions 

(2.2)  (hD+)P  wQn  * 0 , p > 1,  wNn  = g 

are  given.  Then  the  approximation  (2.1),  (2.2)  is  stable  (see  [2])  and 
in  every  finite  time  interval  the  solution  of  (2.1),  (2.2)  converges 
for  h -►  0 to  the  corresponding  solution  of  (1.2)  which  satisfies  the 
boundary  condition 


(2.3)  v(a,t)  ■ g. 

We  new  study  the  behaviour  of  (2.1),  (2.2)  for  t » with  fixed  k,h. 
We  have 

Theorem  2.1.  Let  k,h  be  fixed.  The  solutions  of  (2.1),  (2.2)  con- 
verge exponentially  fast  to  the  steady  state  solution 

(2.1»)  *v=  g , V»0,1,2,...,H 

as  n ♦ •. 
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Proof.  We  write  (2.1),  (2.2)  ih  matrix  form 

(2.5)  71+1  = Qw1+b,  b = (0,0,. ..,0,8)* 

where  w * (w^, . . i ^w^) ' is  a column  vector  and  Q an  N+l  N+l  matrix. 
It  is  well  known  that  the  solutions  of  (2.1),  (2.2)  converge  to  the 
steady  state  Solution  (2.U)  if  and  only  if  the  eigenvalues  z of  Q 
satisfy  |z|  < li  The  eigensolutions  <f>  have  the  form 


v V 

Vi  * Vz  ’ 


where  tc^,  <2  are  the  roots  of 

(2.6)  ac  = K + ~ \ (k2-1)  + \ X2(k-1)2 


<p  must  satisfy  the  homogeneous  boundary  conditions  (2.2),  i.e. 

‘’l^l"1^  + a2^2_1^P  = °»  aiKlS  + °2K2N  = 0 * 

This  system  has  a nontrivial  solution  if  and  only  if 

(2.7)  (Krl)P  </  = k1N(k2~1)P  . 

An  easy  calculation  shows  that  for  sufficiently  small  h the  relations 

(2.6),  (2.7)  imply 

|z|  < 1 - 6,  5 > 0 independent  of  h 

for  all  eigenvalues  of  Q.  This  proves  the  theorem. 

We  next  investigate  the  case  where  we  use  extrapolation  also  at 
x ■ a,  i.e.  we  replace  (2.2)  by 

(2.8)  hD  v.D  ■ hD  w_  ,n  ■ 0 . (For  simplicity  we  set  p*l.) 

In  this  case  we  have 
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Theorem  2.2.  Assume  that  v(x,0)  * f(x)  is  a smooth  function  vith 
df(a)/dx  = 0.  Then  the  solutions  of  (2.1),  (2.8)  converge  for  h -*■  0 
to  the  solution  of  (1.2)  which  satisfies  the  boundary  condition 

(2.9)  v(a,t)  “ f(a) 


Proof.  Let 


c 


c 


Vvn  a D+Vvn,  V=0,1,2,...,N-1. 


Then  v^n  is  the  solution  of  (2.1)  which  satisfies  the  initial  and 
boundary  conditions 

vv°  * D+fv  , V=l,2,...,N-2;  wQn  = w^11  = 0,  n=0,l,2,.., 

S.  Parter  [3]  has  shown  that  the  w^n  are  uniformly  hounded  and  for 
any  6 > 0 


C 


.0 


(2.10)  sup  |y(x..,t  )-w,n|  0 as  h -*■  0 . 

6<x^< a , 0 <t T v n v 

Here  y is  the  solution  of  (1.2)  with  initial  and  boundary  conditions 
y(x,0)  = df(x)/dx  , y(a,t)  = 0 * df(a)/dx  . 

Therefore  we  have  also  that 


lim  D-W. 


b+0 


CTN-1  = y(a«tn)  = 0 • fim  hD+D-WH-l  " 0 

h-*0 


and  by  (2.1) 


low 


n+1  _ n 

lim  — r * 0 , i.e.  lim  w n ■ f(a) 

h+0  * fr+0  N_1 


v-j 


vj 


and  by  (2.10) 
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sup  |v(x  ,t)-v(x  ,t)|  -*•  0 as  h -*■  0, 

(2.11) 

a 

v(x,t)  « f (a)  - / y(£,t)d£;  . 


v(x,t)  is  the  solution  of  (2.1)  which  satisfies  the  boundary  condition 
(2.9).  Also,  for  0 <_  x^  <_  6 ve  have 


(2.12) 


|v(xv,t)-w(xv>t)|  < | v(x,t)-v(x,t)|  + |v(xv,t)-v(x,t)| 
+ |w(x,t)-w(xv,t)|  < const.  5 , 


where  x = x^  is  a gridpoint  with  <5  <_  £ 5 + h. 

theorem. 

For  steady  state  calculations  we  have 


This  proves  the 


Theorem  2.?.  For  fixed  k,h  the  scheme  (2.1),  (2.8)  is  weakly 

00 

convergent  as  n -*■  <*>  i.e.  the  limit  function  w depends  on  the 
initial  values 

Proof.  We  write  (2.1),  (2.8)  again  in  the  matrix  form  (2.5).  An 
easy  calculation  shows  that  z **  1 is  a simple  eigenvalue  corresponding 
to  the  eigenfunction  = const.. All  the  remaining  eigenvalues  satisfy 
again  | z|  <_  1 — 6 < 1.  This  proves  the  theorem. 

Thus  the  steady  state  depends  on  the  initial  function  f(x).  If  we 
replace  f(x)  by  f(x)  + const,  then  also  the  steady  state  changes 
by  this  constant.  Observe  that  the  rate  of  convergence  is  the  same  as 
in  theorem  2.1. 

As  a third  alternative  ve  consider 


(2.13) 


“♦Vi 


o , 


i.e.  the  boundary  values  are  given  on  the  "wrong  side".  This  corresponds, 
for  example,  to  the  subsonic  case  in  fluid  mechanics  where  all  the 
variables  are  prescribed  on  the  inflow  side.  In  a similar  way  as  theorem 
2.2  on  can  prove 
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Theorem  2.4.  Assume  that  the  conditions  of  theorem  2.2  are  satis- 
fied. Then  the  solution  of  (2.1),  (2.13)  converges  for  h -*•  0 on  any 
finite  domain  6<_x£a,  0£t£T  to  the  solution  of  (1.2),  (2.9). 

For  steady  state  Calculations ve  have 

Theorem  2.5.  For  fixed  k,h  the  solution  of  (2*1),  (2.13)  converges 

00  4 

for  n -*■  00  to  the  unique  steady  state  w^  = g,  v=0,l,2,. W ,N  if  and 
only  if  N is  odd.  Also,  the  speed  of  convergence  is  extremely  slow. 

Proof.  The  condition  (2.?)  for  the  existence  of  a nontrivial 
solution  becomes  now 

(2.14)  (^-1)*/  1 = k2N-1(k2-1)  . 

Using  (2.6)  and  (2.14)  an  easy  calculation  shows  that  there  is  one 
eigenvalue  of  the  form 

Z S 1 + (‘1)NX  M (I^)N_1  • X=k/h* 

while  all  the  other  eigenvalues  satisfy  |z|  <_!  — 6 < 1.  Therefore,  we 
obtain  growing  solutions  if  N is  even,  but  converging  solutions  if 
N is  odd.  The  rate  of  convergence  in  the  last  case  is  determined  by 
the  above  eigenvalue,  i.e.  it  is  extremely  slow.  Furthermore,  the 
solution  will  contain  oscillatory  components  near  the  boundary. 

For  converging  solutions  the  steady  state  is  given  by 

(2.15)  wv  = a1<1  + o2<2 
where  K1»IC2  811:6  the  roo1;s 

K2  - 1 + X(k-1)2  = 0,  i.e.  K±  = (l-A)/(l+A),  <2*1 

vv  must  satisfy  the  boundary  conditions  (2.13).  Therefore 

°1  + °2  = 6*  * 02k2N~1(k2-1)  * 0 

which  gives  us  = 0,  a2  * g and  the  theorem  follows. 


By  theorem  2.h  the  solution  of  (2.1),  (2.13)  converges  to  rr • 
solution  of  (1.2),  (2.9)  on  every  finite  time  interval.  Since 
v(x,t)  = f(a)  for  t a this  seems  to  indicate  a contradict ir: 
the  result  of  theorem  2.5*  However,  it  can  be  explained  in  the 
following  way.  For  fixed  values  of  k,h,  a typical  calculation  - 
after  a relatively  short  time  give  a good  approximation  of  v(x,t 
f(a),  but  it  will  later  change  very  slowly  and  finally  converge  t 
The  following  calculation  illustrate  this  point. -We  compute  th» 
of  (2.1),  (2.13)  on  the  interval  0 £ x £ 1 using 


fv  = sin  xv 


9 


hD+WN-l 


0 


as  initial  and  boundary  conditions.  The  figures  show  v(2/3,t ) as  . 
function  of  time  in  two  different  scales.  It  is  very  tempting  to  ex 
sider  the  solution  as  a steady  state  after  a relatively  short  ti*.' 
However,  as  can  be  seen  from  fig.  1,  the  solution  at  t * 2 has  notr  . 
to  do  with  the  final  steady  state. 

We  can  also  use  higher  order  extrapolation  at  the  boundary 

(2.16)  (hD+)P  wQn  = (hD+)P  wN_in  = 0 , p > 1. 


The  condition  (2.7)  for  a nontrivial  solution  becomes 
(x1-1)p(x2-1)p(k2N_p-x1N"p)  = 0 


and  there  is  an  eigenvalue  z = 1 of  order  p which  has  only  or.«  * r 

. p-1 

solution.  Therefore  there  are  solutions  which  grow  like  t and  for 
steady  state  calculations  no  convergence  will  occur. 

For  calculations  on  a finite  time  interval,  we  can,  using  the 
technique  as  in  theorem  2.2,  prove  that  the  solutions  converge  for 
h -*■  0 to  a solution  of  (1.2)  which  satisfies  the  boundary  condition 

3Pv(a,t)/3xP  * 0 . 

The  results  in  this  section  can  also  be  generalized  to  equations 


3v/3t  « a 3v/3x 


a > 0 


end  difference  schemes 

vvn+1  = (l-wkD0+Bk2D+D_)wvn  . 

Stability  for  the  Cauchy  problem  requires 

0 <_  —■  , X2  <_  20  where  X = ak/h. 

We  note  in  particular,  that  for  the  case  of  overspecification  at 
x * 0,  we  get  convergence  to  a steady  state  if  and  only  if  at  least 
one  of  the  following  conditions  is  satisfied. 

1)  N is  odd,  2)  X < 20  . 

Therefore,  by  strengthening  the  von  Ifeumann condition,  we  can  avoid  the 
restriction  on  N.  However,  the  speed  of  convergence  is  still  extremely 
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3.  Difference  approximations  for  systems  with  constant  coefficients. 
In  this  section  ve  will  make  a few  comments  on  the  numerical  solution 
of  systems 


3v/3t  = A 3v/3x 
r~(D  ~(n)\. 


, 0 <.  x <_  a 


where  v * (v'4-7.,, . . ,vVIi7)'  is  a vector  function  with  n components 
and  A a constant  n*n  matrix.  We  assume  that  A can  he  transformed 
to  diagonal  form  and  that  the  eigenvalues  of  A are  real  and  nonzero, 
i.e.  there  exist  a nonsingular  transformation  S such  that 


(3.1) 


where 


A ■ S-1A  S * 


(3.2) 


*1 


f\  0 


0 X2  0 . . . 0 


fX  r+1  0 • • • • 


0 X 


> 0,  A, 


r+2 


0 ‘ * 0 < 0 


V 


0 X 


c 


i o 


are  positive  and  negative  definite  diagonal  matrices  respectively. 
The  Lax-Wendroff  scheme  is  given  by 


(3.3) 


~ n+1 

Vv 


V - (l+kAD0+|k2A2D+D_)wvn  , 0 < V < N 


To  determine  its  solution  we  have  again  to  specify  boundary  conditions 

at  x * 0 and  x ■ a.  If  the  boundary  conditions  are  of  the  same  type 

«*• 

for  all  components  of  w,  for  example 


(3.U) 


wQn  « g , hD+wN_1D  « 0 , 


then  all  the  results  of  the  last  section  apply.  We  need  only  to  introduce 

,-l~ 


new  variables  • S vv  to  obtain 


l 

I 
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L 


O 

4 

. o 
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wvn+1  «=  (l+kAD0+|k2A2D+D_)wvn  , 0 < V < N 

(3.5) 

□ . _ n 

wQ  ■ g,  hD+wR_1  = o 

i.e.  n decoupled  problems  of  type  (2;l),  (2.13). 

This  transformation  can  also  be  used  to  discuss  general  boundary 
conditions.  For  example,  specify  at  x = 0 correct  boundary  conditions 
and  extrapolate  at  x = a all  variables.  Then,  after  transformation, 
these  boundary  conditions  become 


<»on>"  * ao(V>°  ♦ • (“’.)P(''oI)°  “ ato.l'H11' 


Ln  , , II  «n 


I*n 


\P 


(3.6) 


“+Vi  = 0 • 


Here  v*,  v11  correspond  to  the  partition  of  A and  R,  Q denote 
(n-i&r  and  r^n~r)  matrices  respectively.  Then  v converges  to  a 
solution  of 


(3.7)  3v/3t  * A 3v/3x 
with  boundary  conditions 

(3.8)  vn(0,t)  * RgV^O.t)  + gon(t). 

By  (3.7)  the  relation  3v*/3x  * 0 implies 


3v*(a,t)/dx  * 0 . 

3v*/3t  * 0 and  therefore 


(3.9)  vI(a,t)  = vI(a,0)  . 

Thus  in  steady  state  calculations  the  steady  state  depend  on  the 
initial  values. 


• 1 
li 


( 


c 


c 


I . • 

■ O 

I • 


r 

j 

» J 
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U.  General  considerations  for  problems  in  one  space  dimension.  In 
this  section  we  consider  hyperbolic  systems 

3u/3t  = A(x,t)3u/3x 

with  variable  coefficients  in  the  quarter  space  0 £ x < « , t > 0 
where  the  matrix  A depends  smoothly  on  x,t.  Ve  assume  that  there  is 
a smooth  transformation  S(x,t)  such  that  (3.1),. (3.2)  hold  for  every 
fixed  x,t.  Then  we  can  introduce  new  dependent  variables  u = S \i  and 
obtain 

(U.l)  3u/3t  = A 3u/3x  + Bu,  B * -S_13S/3t  + S_1A  3S/3x  . 

For  t = 0 initial  values 

(U.2)  u(x,0)  * f(x)  , 0 £ x < ®> 

and  for  x * 0 boundary  conditions 

(1».3)  uII(0,t)  = RgU^O.t)  + g0i:E(t)  , t > 0 

are  given.  Here  f(x),  g^^f)  are  smooth  functions,  u*  ■ (u^ , . . . ,u^ ) ' , 

U11  * (U(r+1),...,u(n))'  correspond  to  the  partition  of  A , and  RQ 

is  an  (n-r)xr  matrix.  Thus  the  number  of  boundary  conditions  is  equal 
to  the  number  of  characteristics  with  enter  the  region  x > 0,t  > 0. 

If  we  want  to  solve  the  above  problem  numerically  we  have  to  replace 
the  infinite  interval  0 £ x < °»  by  a finite  one.  There  are  two  ways 
to  do  this. 

1)  Determine  some  transformation  of  the  independent  variable  x 
which  transforms  0 £ x < • into  0 £ x £ a. 

2)  Replace  0 £ x < ® directly  by  0 £ x £ a and  find  at  x * a 
a function  g^I(t)  such  that  the  solution  of 

(k.U)  3v/3t  ■ A 3v/3x  + Bv,  0£x£a,  t > 0 

with  initial  conditions 

(fc.5)  w(x,0)  » f(x)  , 0 £ x £ a 

. 

• . 


- - - - 
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and  boundary  conditions 

(4.6)  vII(0,t)  - RQvI(0,t)  ♦ g0IX(t)  , vI(a,t)  « g^t)  , t > 0 

differs  at  most  slightly  from  the  solution  of  the  original  problem 
(4.1)— (4 .3) • 

Connected  with  the  above  equations  is  the  following  problem. 
Consider  the  system  (4.1) 

(4.7)  3y/3t  »=  A 3y/3x  + By,  x>  a,  t >,  0 
for  x > a with  initial  values 


L 


(4.8) 


y(x,0)  * f(x)  , x > a 


and  boundary  conditions 

(4.9) 


II,  ..  II,.  , 
y (a,t)  = gx  (t)  . 


We  can  express  the  solution  u(x,t)  of  (4.l)-(4.3)  with  help  of 
v(x,t)  and  y(x,t).  The  following  lerana  is  obvious. 

Lemna  4.1. 


. c 


Q 


(4.10) 


u(x,t)  ■ J 


v(x,t)  , X < a 
y(x,t)  , x > a 


» t 0, 


if  and  only  if 

(4.11)  vI(a,t)  ■ g1I(t)  ■ yI(a,t),  vXI(a,t)  * g1II(t)  *=  yIJ(a,t)  . 

Ve  can  state  this  result  in  another  way.  Let  f(x)  be  fixed.  Then 

the  problem  (4.7)-(4.9)  has  for  every  g.II(t)  a unique  solution 

X ^ XI 

y(x,t).  In  particular,  y (a,t)  is  uniquely  determined  by  y (a,t)  ■ 

gj*Z(t).  Thus  there  is  a linear  operator  R^(t ) such  that 

(4.12)  yT(aft)  - K1(«)yIZ(att)  ♦ b*(t)  . 
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c 

c 

• C 

t * * 

! P 


Here  b*(t)  is  determined  by  f(x)  and  H_ (t ) depends  on  y (a, 5) 

I 

for  0 £ £ £ t.  Assume  that  R1  and  b are  known.  Choose  the  function 
g^(t)  in  (4.6)  such  that 

(**.13)  vI(a,t)  * R^tjv^ta.t)  + bI(t)  . 

We  have 

Lemma  4.2.  u(x,t)  * v(x,t)  for  0 <_x.<  a,  t ^ 0 if  and  only  if 
(4.13)  holds.  Thus  (4.13)  can  be  considered  as  the  missing  boundary 
condition. 

Proof.  Assume  that  v(x,t)  satisfies  (4.13).  Choose  g^**(t)  in 
(4.9)  such  that  y**(a,t)  = g^*(t)  * vII(a,t).  Then  by  (4.12)  and  (U.13) 

also  y*(a,t)  = v*(a,t)  and  by  lemma  U.l  we  obtain  u(x,t)  = v(x,t) 
for  0<  x < a,t  0.  Conversely,  if  u(x,t)  is  a solution  of 
(U.l)-(U.3)  then  it  is  also  a solution  of  (U.7)-(U.9)  and  must  satisfy 
(U.12).  Therefore,  if  u(x,t)  * v(x,t)  for  0 £ x£  a,t  5^  0 then 
v(x,t)  must  staisfy  (U.13).  This  proves  the  lemma. 

In  general  R.(t),b*(t)  are  very  complicated.  There  are  only  some 

j 

special  cases  where  R^  and  b can  be  represented  in  a simple  way. 

1)  r = 0,  i.e.  all  eigenvalues  of  A(a,t)  are  negative,  i.e.  all 
characteristics  at  x * a point  out  of  the  region  0 _<  x _<  a,t  >.  0.  In 
this  case  y**  * y and  the  relation  (U.12)  is  empty.  Therefore  v(x,t) 
does  not  need  to  satisfy  any  boundary  conditions  at  x = a. 

2)  r > 0 but  A(x,t)  = A(a)  is  a constant  matrix  for  x >_  a.  Then 
the  transformation  S is  independent  of  x,t  and  B = 0 for  x >_  a,. 
Therefore  we  can  compute  the  first  r components  of  y(x,t)  explicitely. 
They  are  given  by 

y^(x,t)  ■ f^(x+A.t)  , j»l,2,...,r; 

J 

leading  to 

y^(a,t)  ■ f^(a+A.t)  , j-1,2,. . . ,r. 

J 


J 
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r 

c 


o 
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The  relation  (U.12)  holds  with  R-  **  0 and  b*  * (f^(a+X-t ) , . . . ,f ^(a+X  t))’.; 

l 1 r ; 

It  is  clear  that  v(x,t)  depends  very  much  on  the  initial  values 
for  x ^ a.  There  is  only  one  simple  case,  namely 

f(x)  = f = const,  for  x > a . 

09 


Then  ( U .13)  becomes 

(U.lU)  vI(a,t)  * f1^  .m) 


(U.lU)  is  equivalent  with 


(U.lU)'  3v*(a,t)/3x  * 0 , 


because  ( U . lU ) * implies  3v*(a,t)/3t  *>  0,  i.e.  v*(a,t)  = vI(a,0)  = f^  . 
Thus  the  boundary  conditions  (3.**)  (extrapolation  of  all  the  variables 
at  x = a)  can  be  used. 

Without  restriction  we  can  assume  that  - f =0.  Otherwise  we  would 

00 

• . ~ I • 

consider  the  function  u a u - f . Then  also  b * 0 and  the  relation 

09 

(U.13)  becomes 

vI(a,t)  ■ 0 , for  t ^ 0 . 


The  last  relation  means  that  we  set  the  characteristic  variables 
associated  with  the  "ingoing”  characteristics  equal  to  zero,  or  as 
B.  Engquist  and  A.  Majda  [1]  call  it,  that  no  reflection  takes  place 
at  the  boundary.  Therefore  this  principle  is  useful  if  we  subtract  from 
the  solution  its  constant  state  at  infinity. 

From  now  on  we  shall  always  make 

Assumption  U.l.  f(x)  = 0 for  x >_  a,  i.e.  b*(t)  * 0. 


Assume  that  A depends  on  x,t  also  for  x >_  a.  Then  B ^ 0 and 
in  general  y*,  y**  are  coupled  and  we  cannot  determine  R^(t)  without 
making  detailed  calculations  of  y(x,t).  Only  if  A(x,t)  converges  to 
a constant  matrix  AQ  as  x -*•  « we  can  do  better.  We  make 


iT 


▼I(a,t)  can  also  be  obtained  in  a fairly  simple  way  if  we  require 
only  that  B has  the  form 


i 


» 

I 


i 
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i 
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Assumption  U.2.  The  matrix  B(x,t)  can  he  written  in  the  form 
(U.15)  B(x,t)  * ({((xjBj^tx.t) 

where  $(x)  is  a scalar  function  with 

(U.16)  / |$(x)| dx  <_  c 

a 

and 

(U.17)  | |b| | * sup  |B(x,t)|<l. 

a<x<»,t>0 

Here  | B^(x,t ) | denotes  the  maximum  norm  at  the  point  x,ti 
Assume  for  example  that  A is  of  the  form 

A(x)  = Aq  + — A.  (x,t)  , x>a 

x 

and  that  the  eigenvalues  of  AQ  are  all  distinct.  Then  S is  of  the 
form 

S ■ S + ~ S (x,t)  , i.e.  B « B . 
x ‘ x 

Thus 

00 

4>  * const./x^  , / |<|>|dx  <_  const. /a  . 

a 

We  need 

Lemma  U.3.  Consider  the  system 

3y/3t  * A 3y/3x  + $(x)G(x,t) 

for  x >_  a,  t > 0 with  zero  initial  and  boundary  conditions 
y(x,0)  ■ 0 for  x > a,  yJI(a,t)  * 0 for  t >_  0 . 

Assume  there  is  a constant  Xq  > 0 such  that 

(U.18)  inf  |X.|>A_. 

t,x,j  J 0 


_ 


U 
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Tben 

||y||  < c/A  ||G||  , ||y||  - «up  |y(l)|  . 

x.t.i 


c 

c 


! 


t 

I I 


Proof.  We  need  to  prove  the  estimate  only  for  scalar  equations 
3y/3t  * A(x,t)  3y/3x  + <J>(x)G(x,t) 

vith  appropriate  boundary  conditions.  Using  the  method  of  characteristics 
ve  obtain 


dy/ds  * $(x)G(x,t)  , dt/ds  * 1 , dx/ds  = -A(x,t)  . 


Therefore 

80  80 

|y(x(s  ),t(s  ))|  * | / <f>(x(s))G(x(s),t(s))ds|  < ||G||  / |g(x(s))|ds< 

0 0 

illGll  / |g(x)| |A_1(x)|dx  < (e/A  ) ||G||  , 

0 u 

vhich  proves  the  lemma. 

The  last  lemma  gives  us 

Theorem  U.l.  Assume  that  the  assumptions  U.l  and  U.2  are  valid  and 
that  (U.l8)  holds.  If  c/A^  < 1 then  the  solution  of  (U.T)  can  be  obtaind 
by  the  iteration  proeess 

3yn+1/3t  - A 3yn+1/3x  * 4>(x)B1(x,t)yn 

(U.19) 

yn+l(x»0)  “ 0 » » 

n*0»l»...  ; yQ(x,t)  * 0 . 

froot.  Let  yQ+1  ■ yQ+1  - yQ.  By  lemma  U.3 

I lyn+il  I il<cA0)  llynll  . n»l,2,... 


J 


/ 

i 


I 


and  the  convergence  follows.  This  proves  the  theorem. 


F 
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In  the  usual  way  we  obtain  the  estimate 

-S  — 

1 ly'ri1 1 1 Ji  1 lJl  1 1 *^*0  1 15,2-5,11 1 


<S«cA0» 


and  by  (U.19) 


(U.20) 


y^(x,t)  = 0 

3yXI/at  - a2  3yXI/at  * o , yXI(x,o)  = o,  yXI(a,t)  « g^(t) 


Thus  if  we  allow  an  error  of  order  $(c/Aq)  extrapolation  or  the 
principle  of  no  reflection  at  x = a is  appropriate.  If  this  error 
is  not  acceptable  then  we  have  to  compute  yX  defined  by 

(k.2l)  ay*/  t - Ax  3y^/3x  = ♦ B^yJ1,  y*(x,0)  = 0, 

where  is  defined  by 


Then 


and  we  assume  that  this  error  is  tolerable. 

If  one  is  interested  in  the  whole  time  dependent  process  then  not 
much  is  gained  because  to  compute  y*(a,t)  as  a function  of  g**(t) 
we  have  to  solve  (U.20)  and  (U.21)  completely.  Thus  one  should  instead 
■ake  a so  large  that  an  error  of  order  (J/ic/X^)  can  be  tolerated. 

There  is  one  exception.  If  c|(B12||  is  small,  i.e.  the  in-  and 
outgoing  waves  are  almost  decoupled,  then  y*  is  almost  zero  and  the 
complete  solution  is  essentially  given  by  (U.20).  In  this  case  the 
principle  of  no  reflection  is  appropriate. 

. The  situation  becomes  more  favourable  if  we  are  only  interested 

in  the  steady  state  solution,  i.e.  consider  the  case  that 

' - 

■ g^X(»)  and  B^x.t)  ♦ B^x),  $(x,t)  •*  A(x)  . 

* •-V'— • 
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Then  for  t -*■  ® 


Thus  y2  satisfies  in  the  limit  the  relation 

y^(a,t)  = C g l1  (.<»),  C = / <p  A^1  B^d?  . 

00 


Therefore,  by  lemma  4.2,  we  should  use  for  v(x,t)  at  x = a the 
boundary  conditions 


vI(a,t)  = C vII(a,t)  . 


Observe,  however,  that  one  has  to  know  the  asymptotic  expansion  of  A 
in  detail  to  compute  C. 
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5.  Problems  in  tvo  space  dimensions.  We  start  with  am  example.  Con- 
sider the  system 


(5.1) 


3w 

at 


r+i  o . 

3w 

0 -1/  3x 


0 1 


3w 

1 0/  ay 


on  the  domain  ft  am  given  in  fig.  1 or  fig.  2.  For  t = 0 initial 
values 


(5.2) 


u(x,y,0)  = v(x,y ,0)  , (x,y)£  ft  , t *=  0 , 


C 


end  for  (x,y)^3  ft  boundary  conditions  are  given.  In  particular  we 


assume  that  on  B, 


(5.3) 


v(0»y»t)  = h(y)  , t ^ 0,  x = 0, 


and  on  B^B^ 


(5.1*) 


u(x,y,t)  * 0,  (x,y)  £.  B3  U t >.  0 


C 


O 


Without  restriction  we  can  eissume  that  B^jB^  are  given  by  x a,y  * 1 
and  x > a, y = 0 respectively. 

We  wemt  to  determine  the  solution  of  the  above  problem  on  ft^  only 
Therefore  we  need  one  relation  between  u,v  on  B^.  The  boundary 
conditions  (5.1*)  imply  that  for  x >.  a the  solution  of  our  problem  can 
be  expanded  in  Fourier  series 

00  00 
(5.5)  u ■ l u(x,(i),t)sinrrwy,  v * v (x,t)  + l v(x,a>,t)cosirwy  . 


u»l 


w»l 


Introducing  (5.5)  into  (5.1)  gives  us  for  x £ a 


(5.6) 


3vQ/3t  » -3v0/3x, 
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I 

c 

I * 

C 

4 

A • 
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For  every  fixed  frequency  u>  f 0 the  system  (5-7)  is  of  the  same  form 
as  the  onedimensional  problem  (U.l)  and  the  results  of  the  last  section 
apply.  Therefore  if  ve  are  interested  in  the  whole  time  dependent  process 
not  much  is  gained  by  restricting  the  problem  to  the  domain  except  in 
the  case  that  v(u>,x,t)  can  be  neglected  for  a)=l,2, . . . ;.  Then  the 
solution  of  our  problem  is  essentially  independent  of  y. 

The  situation  is  again  different  if  we  are  only  interested  in  the 
steady  state.  Then  we  can  replace  (5.6),  (5.7)  by 

(5.8)  3vq/3x  = 0 


The  solution  of  (5.9)  is  given  by 


u(x^.t)  - -e^n(x-a);< 


which  gives  us  for  x * a 

(5.10)  u(a,u>,t)  ■ -v(a,w,t)  . 

Therefore  we  have  determined  the  desired  relation 
u(a,y,t)  • Rv(a,y,t) 

because  if  ve  know  v then  ve  can  expand  it  in  a cosine  series  and  by 
(5.10)  we  obtain  the  sine  Beries  for  u.  This  process  can  be  performed 
■Meric ally  by  using  the  fast  Fourier  transform. 

Zt  is  not  necessary  that  the  coefficients  of  (5*1)  are  constant.  Ve 
cam  also  consider  systems 
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3w 

at 


3w 

3y 


where  the  coefficients  depend  on  y.  The  steady  state  solution  can 
again  he  solved  in  terms  of  eigenfunctions 

w * eKX  <j)(y)  , Real*  < 0, 

which  gives  us  the  desired  relation  between  u and  v on  Bg. 

Whether  this  procedure  is  feasible  numerically  depends  on  how  easy  it 
is  to  compute  the  eigenfunctions  and  how  many  are  needed  to  represent 
the  relation  between  u and  v. 

All  the  results  can  be  carried  over  to  general  hyperbolic  systems 
3w/at  = A^3w/3x  + Ag  3w/3y  . 

Again,  if  one  is  interested  in  the  complete  time  dependent  process  then 
not  much  is  gained  by  restricting  the  problem  to  the  domain  ft,  except 
in  the  case  where  we  can  neglect  Ag  3w/9y.  If  we  are  only  interested 
in  the  steady  state  then  we  can  obtain  the  desired  relations  on  Bg 
between  the  components  of  w by  solving  the  steady  state  equations 

A1  3w/3x  + Ag  3w/3y  « 0 

in  terms  of  eigenfunction  expansions. 

We  want  to  point  out  that  there  are  problems  in  two  space  dimensions 
which  can  be  transformed  on  one-dimensional  form.  Consider  for  example 
equation  (5.1)  on  some  domain  ft  containing  the  origin.  We  introduce 
polar  coordinates  by 

x * r cos0 
y * r sin0 

and  obtain  from  (5.1) 

$v  . A0"8  8in9  \ 3.  4 1 
Jt  \.io0  -00,6/  9r  r 


\ 

| 

i 

! 





I 


-5.U- 


A 

I 


* 


Defining  new  dependent  varieties  by 

0 . 9 


cos  2 
6 


sin  g 
6 


u = T(0)u 


sin  g -cos  - 


we  get  the  equivalent  system 


Therefore,  if  the  initial  data  are  such  that  u(r,0,O)  is  independent  of 
0 , then  u(r,9,t)  is  independent  of  0 for  all  t,  t > 0.  In  particular, 
we  get  a nonreflecting  boundary  condition  if  u*  * 0 is  specified  at 
all  points  on  3ft  . 


1 e 

i • • 
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